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Abstract 

An efficient method to describe the nonhnear evolution of Stimulated Brillouin Scattering in long 
scale-length plasmas is presented. The method is based on a decomposition of the hydrodynamics 
variables in long- and short-wavelength components. It makes it possible to describe the self- 
consistent coupling between the plasma hydrodynamics, Stimulated Brillouin Scattering, and the 
generation of harmonics of the excited ion acoustic wave (lAW). This description is benchmarked 
numerically and proves to be reliable even in the case of an undamped ion acoustic wave. The 
momentum transferred from the electromagnetic waves to the plasma ions is found to induce a 
plasma flow which modifies the resonant three wave coupling between the lAW and the light waves. 
A novel picture of SBS arises, in which both JAW harmonics and flow modification reduce the 
coherence of SBS by inducing local defects in the density and velocity profiles. The spatial domains 
of Stimulated Brillouin activity are separated by these defects and are consequently uncorrelated, 
resulting in a broad and structured spectrum of the scattered light and in a temporally chaotic 
reflectivity. 
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The description of parametric instabilities in laser-produced plasmas using simple coupled 
mode equations for three wave interaction is no longer sufficient whenever the longitudinal 
plasma waves are driven to large amplitudes. Then the nonlinearities of the longitudinal 
wave can induce detuning with respect to the three wave resonance. This is one of the 
reasons usually invoked to explain why these simplified models overestimate the scattering 
levels of Stimulated Brillouin Scattering (SBS). In this article we concentrate on SBS, which 
is the process by which the incident laser wave couples to an ion acoustic wave (lAW) to 
give rise to a scattered transverse wave. The generation of the harmonics due to the lAW 
fluid-type nonlinearity is already known to be able to reduce significantly 

the SBS refiectivity when compared with the results involviiig_simply a linearized lAW. 
However, the previous fiuid-type models for SBS in Refs. UQ, Q , aimed at taking into 
account the lAW nonlinearity, were incomplete because they did not properly describe the 
fiow modification j^, Q] caused^ by the incident transverse wave momentum deposition. All 
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the mentioned models [l|, 12, 121 14, l3] also ignored multidimensional effects. On the other 
hand, kinetic effects associated with particle trapping [8] give also rise to a nonlinear lAW 
frequency shift and therefore modify the SBS nonlinear behavior. 

In the present Letter, we reconsider the effect of the lAW nonlinearities on SBS by ac- 
counting properly for the flow modiflcation caused by SBS. We flrst derive approximate 
equations describing simultaneously the plasma hydrodynamics (i.e. the long wavelength 
density and flow proflles), SBS, and the harmonic generation of the excited lAW resulting 
from fluid-type nonlinearity. Our method consists in decomposing the fluid variables into 
long and short wavelength components, the latter corresponding to the SBS generated lAW 
and its harmonics Our new code, based on this harmonic decomposition method, makes 
it possible to describe plasmas of spatial sizes of the order of realistic laser produced plasmas 
(of mm-size, typically), because it does not resolve the lAW /im-scale. We then continued 
a step further by checking the capacity of our approach to account for kinetic effects ef- 
fects by implementing in the lAW propagator a nonlinear frequency shift modeling particle 
trapping. [S| 

The transverse electric fleld is described by E{-x,t) = e"""^"* (^i?+e*^°^ + E_e~^^°''^ + c.c. 
where i?+(x, t) and -E'_(x, t) are the forward- and backward propagating light fleld compo- 
nents, respectively, both enveloped in time and space with respect to the light frequency uq 
and the wave number /cq- This wave number is taken for a flxed reference plasma density 
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Neq which yields, using the critical electron density nc, = to'o(l — Neq/ric) / . For the 
plasma density n{'x,t) and the velocity v(x, t) we use a decomposition separating the long- 
wavelength components A^o(x;^) and vo(x, t) and the short-wavelength components np(x, t) 
and Vp(x, t), with \p\ = 1, 2, 

n = No+ (nie*'""^ + n2e''''' + .. + c.c.) , 
V = vo + (vie*''^^ + vse'^^^ + .. + c.c.) , 

the first (p=0) representing the hydrodynamic evolution, and the terms with p >0, the 
fundamental ion acoustic wave, p =1, excited by SBS, and its harmonics, p >1. The 
reference wave number for the lAW is the wavenumber of backscattering, kg = 2ko, for 
which the ponderomotive force is proportional to oc E^E'^exp{i2koz). 

We use the paraxial approximation to reduce the wave equation for the total electromag- 
netic field E to two "paraxial" equations for i?+(x, t) and i?_(x, t), 

CpUE+) = -i{uJo/cN,^) [niE^ +(iVo - K^)E+] , (1) 
C-^.r{E^) = -t{uJo/cN,^) [nlE+ +{No - N,^)E_] , (2) 

with the paraxial operator £par(-E±) = [dt + c±dz + vt —i{,c^ /2ujq)V']]E±, where c+ and c_ 
stand for the group velocity of the forward/backward propagating light, respectively, with 
c+ = c^ko/ujQ = — c_, and Ut denotes the damping of the transversal waves. The right-hand- 
side (rhs) source terms in equations (0) and (j2)) account for (i) resonant 3-wave coupling due 
to SBS, with the fundamental ion sound wave, rii, and for (ii) refraction on long- wavelength 
density modifications, A^o ~ ^eq, causing e.g. self-focusing. In comparison with the full wave 
equation without decomposition into E±, this model allows a considerably coarser spatial 
resolution and thus much less numerical expense. 

For the long-wavelength hydrodynamic component we use the following set of equations, 
assuming isothermal conditions, and written in the conservative form on the left-hand side 
(Ihs): 

dtNo + ViVoVo = {dtn\^^ , (3) 
dt (iVoVo) + V (iVoVoVo) + clVNo = (4) 
-NoclVUo + (9tnv)jAw , 

where the rhs source terms, {dtn)iAw and (9tnv)iAw, describe the momentum transfer into 
the flow due to the lAW excitation by SBS, with (i9tnv)j^^ = 2cs {2h'si — vq ■ V) \ni\'^/NQ, 
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and (c^i^)iAw = ~2csV (|nip/A^o) • The ponderomotive force is given by VUq = eoV(|-E'+p + 
|i?_P)/r2cTe. The equations describing the lAW driven by SBS, rii, and its harmonics, riiyi 
(using the convention n_/ = for the complex conjugate) can be written, in the so-called 
weak coupling regime, as follows 



with Qi = {l/2)Y.nhni^h for /i 7^ and / ^ h, where c, = [{ZTe + 3Ti)/Mi]i/2 jg ^j^g j^w 
speed (with Z and Mi as the ion charge and mass), vqz the ^-component of the flow vq; Vgi 
and uJi denote the group velocity and the "local" frequency of the l-th lAW harmonic, both 
accounting for the dispersion due to Debye shielding increasing with the harmonic order. 
They are given by Vgi = Cs(l + t^klX^)'^/"^ and by loi{z) = uJs{lksCs) + lksVoz{z) with the 
lAW frequency iOs{k) = kcs{l + k'^Xjj)^^^'^. Equations (P)-® describe what we call the 
harmonic decomposition model. They form a closed system describing SBS in a temporally 
and spatially evolving plasma. They can be shown to conserve momentum 0] at the lowest 
order in l/{ksi\\) and in l/{ksi±y (with the inhomogeneity length i\\ = \dzVo/vQ\~^ and 
i± = |V^t;o/t;or'). 

In the following we emphasize the particular importance (i) of the SBS-induced flow 
modification, originating from the rhs term of Eq. Q as well as of the term (9jnv)iAw 
on the rhs of Eq. and (ii) of the lAW harmonic generation described by the coupling 
terms oc Z^^/i'^m-Zi in the rhs of Eq. (0). In order to stress the effect of each mechanism, 
we neglect for simplicity the lAW damping, (while being aware that the lAW damping 
coefficient is usually of the order of a few percent of the lAW frequency). Indeed, the 
SBS-induced flow modification due to momentum transfer, first pointed out by Rose in 
Ref.j^, cannot be ignored in the regime of absolute instability corresponding to weak lAW 
damping, because it is just in this regime that the stationary ID limit of Eqs. and Q 
exhibits the most pronounced fiow modification. Namely, the generation of the backscattered 
light gives rise to a transfer of momentum to the bulk plasma in the spatial domain of 
SBS activity. This momentum transfer results in a decrease Av = vo^out — "Wo.m < of 
the fiow Vq in the direction of propagation of the laser, the net fiow decrease being given 
by Av ~ -2i?sBs(2eo|^+|V^egre)(l - N^q/2nc). Here, i?sBS denotes the SBS refiectivity 
corresponding to the considered SBS active region. 



\ + I'si + i^i + {vqz + Vgi)dz - i{cs/2lks)V']_ rii 




(5) 
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FIG. 1: SBS reflectivity -RsBS versus time for the case of an undamped lAW with the parameters 
II = 2.5- lO^^W/cm^ for Aq = 1.064/xm at = IkeV, iVo/^c =0.1 (taken at center), 2koXD =0.27, 
Lini ~ 160Ao. The sohd Une is obtained from the decomposition code considering ah terms, the 
dashed hne from the decomposition code disregarding higher lAW harmonics, and the dash-dotted 
hne from the "full" code. 

We have performed simulations on the basis of equations (0)-© and expanded the JAW 
up to its 3rd harmonic, resulting in a set of equations for rii, n2, and n^, with the rhs terms 
Qi = n2n\ + n^n2*, Q2 = n{ + 2n^ri[, and = 3n2ni. We did not observe any significant 
changes when harmonics above the 3rd order were retained, while restricting to less than 3 
harmonics led to important differences. At this stage of our study we restricted ourselves to 
one-dimensional (ID) simulations in order to benchmark our harmonic decomposition code 
against a "complete" ID code which does not make the decomposition corresponding to 
Eqs. (l)-(3). This latter code solves Helmholtz's equation for the total electric field E{z,t) 
on the first hand, and the system of fluid equations for continuity and momentum, with the 
complete ponderomotive force, V\E{z,t)\'^ , as a source term, on the second hand. Here, in 
ID, the operator V reduces to the partial derivative e^dz- 

To ensure equivalent boundary and initial conditions we have considered a realistic case 
similar to an "exploding foil", where an initially heated plasma expands starting from an 
almost box-like density profile, with smooth shoulders, in the interval zi < z < Z2 along 
the laser axis. The plasma profile, with the initial plateau width Li^i ~ I6OA0, successively 
undergoes rarefaction from each side, so that the velocity profile eventually tends to a 
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monotonous curve varying from negative to positive values with vq = in the center. The 
simulation box is chosen in such a way that the rarefaction of the profile does not significantly 
change the boundary conditions for the light fields at the entrance {zent = ^ < ^i) and the 
rear side (2:rear > Z2). The total box size is ^rear = 2000//co — 320Ao, where Aq = 2n/ko 
denotes the laser wavelength. The boundary condition for the incident light at 2; = is a 
constant, -E+(0) = const, whereas the backscattered light is seeded with a noise source at 
the level {\E_{z = Zrear)P) ~ 10~^\E+{z = 0)p and with a spectral bandwidth sufficiently 
larger than the lAW frequency, in order to cover all possible SBS resonances in the profile. 
In the density profile wings left and right of the central plateau (for times t < Li^i/2cs), the 
plasma is strongly inhomogeneous in velocity and density so that SBS is inhibited by the 
strong fiow gradient. 

We carried out our simulations in the absolute instability regime of SBS with undamped 
lAWs, both to examine the role of fiow due to momentum transfer, and to benchmark the 
robustness of our decomposition code. Notice that in the case of completely undamped 
lAWs the SBS saturation level is, according to Refs. 0, [ill, independent of the noise 
level. For the chosen electron density Neq/ric = 0.1 and for the plasma length indicated 
above, the standard three-wave interaction model for undamped lAWs predicts a 

steep increase in the SBS refiectivities -Rsbs as a function of the laser intensity, varying from 
RsBS << 1 for small laser intensities to -Rsbs — 1 for normalized laser intensities above 

= eolEol"^ / ricTe ~ 0.003, with Eq denoting Eq = E+{z = 0). Our simulations comparing 
the decomposition code and the "complete" code show very good agreement, even for the 
extreme case shown in Fig. 1, corresponding to the plasma parameters mentioned above and 
to Og = 0.025 and 2kQX£) = 0.27 (corresponding to an electron temperature Tg = IkeV and 
a laser intensity II = 2.5 ■ lO^^W/cm^ at Aq = 1.064/im), for which the refiectivity would be 
99% in the absence of any lAW nonlinearity or fiow modifications. For lower intensity values, 
the agreement is even more striking. This excellent agreement between the two codes gives 
us confidence in the robustness of the harmonic decomposition description. In the simulation 
presented in Figs. 1 and 2, the maximum amplitudes of the harmonics remained below the 
validity condition for harmonic expansion, namely \ni/N(,q\ < Q^^'^{lks\£,Y for / = 1,2,3. 

It can be observed, in the spatial profiles shown in Fig. 2 for the backscattered inten- 
sity l-E-P, the fundamental lAW amplitude |ni|, the fiow fo, and the plasma profile Nq, 
that the lAW behavior and fiow modifications are entirely connected with the existence of 
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FIG. 2: Spatial profiles, taken at 2koCst = 220, of the backscattered intensity (upper subplot), 
the fundamental lAW square amplitude the flow velocity vq, and the density A'^o inside the 

shoulders of the exploding foil profile. 

"defects" in these spatial profiles and with a non-monotonous character in space (see also 
Ref. 0]). Namely, SBS develops in distinct spatial domains, interrupted by phase defects, 
which originate in the density profile shoulders corresponding to the low density plasma on 
the laser entrance side, and which then propagate into the profile plateau. Thus the SBS 
activity in each spatial domain appears to be uncorrelated, due to their different origin in 
the inhomogeneous velocity profile vq- This feature reflects in the structured nature of the 
backscattered light temporal spectrum, shown in Fig. 3 in which distinct peaks appear, and, 
consequently, in the temporally chaotic behavior of the reflectivity. 
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Our decomposition description makes it possible to discriminate the relative importance of 
the various effects contributing to the nonstationary behavior in the SBS reflectivity, as seen 
in Fig. 1. By suppressing parts of these effects in different runs, we have found that the most 
important effect is the excitation of the lAW harmonics: namely, retaining the harmonic 
excitation and neglecting the SBS-induced flow modification lead to results that remain in 
reasonably good agreement with the exact model, whereas, retaining the flow modification, 
but ignoring the harmonics leads to unphysically high levels of lAW amplitudes. It follows 
from these observations that a realistic modeling of SBS requires the proper description 
of the lAW harmonics. The ID simulations presented here would correspond to the SBS 
development in a long laser hot spot. We have recently carried out 2D simulations which 
confirm the relevance in 2D of the scenario described above whenever the hot spot focus 
is not far (less than approximately one Rayleigh length) from the transition between the 
inhomogeneous and the homogeneous domain of the plasma density profile (i. e. the shoulder 
of the expanding plasma in our case) . 

Increasing the laser intensity induces stronger lAW amplitudes at which ion and/or elec- 
tron kinetic effects take place. We have included phenomenologically weak ion kinetic 
effects in our decomposition model by adding a non-linear frequency shift of the form 
—iri\ni/Neq\^/'^^ jS] in the propagator appearing in the Ihs of Eq. (0) describing the evo- 
lution of the lAW fundamental component and of its harmonics. Although this is subject of 
work in progress, let us mention that we have solved numerically the corresponding equation 
Eqs. (IH)-©, and we find that for a positive and sufficiently large rj coefficient (~ 0.5 . . . 0.7), 
this shift can smooth out the effect induced by the harmonics and the flow, in a way such that 
(i) the "defects" are less pronounced and (ii) the SBS reflectivity diminishes, but without 
exhibiting a strong nonstationary behavior. In conclusion, we have shown that the SBS mod- 
eling presented here, based on a harmonic decomposition of the hydrodynamics variables, 
represent a promising way to describe laser plasma interaction in long scale-length plasmas. 
We have benchmarked our code based on the harmonic decomposition in the extreme limit 
of the absolute instability regime by neglecting the lAW damping. A novel picture of SBS 
arises in which an incoherent superposition of scattered light generated in distinct spatial 
domains in the velocity profile leads to a nonstationary character of the SBS reflectivity and 
to a significant reduction in the time averaged reflectivity. This harmonic decomposition 
description appears to be sufficiently robust and versatile to allow further sophistication by 
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FIG. 3: Spectrum of the backscattered light, corresponding to Fig. 1. The frequency is shifted 
with respect to the incident laser frequency loq. 

including additional mechanisms such as kinetic effects via an amplitude-dependent nonlin- 
ear frequency shift. We currently work on a generalization of the harmonic decomposition 
method in order to include the subharmonic lAW decay. 

The numerical simulation were carried out thanks to the access to the facilities of IDRIS 
at Orsay, France. The authors would like to acknowledge fruitful discussion with L. Divol, 
J. Myatt, C. Riconda, H. A. Rose, and W. Rozmus. 
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The description of parametric instabilities in laser-produced plasmas using simple coupled 
mode equations for three wave interaction is no longer sufficient whenever the longitudinal 
plasma waves are driven to large amplitudes. Then the nonlinearities of the longitudinal 
wave can induce detuning with respect to the three wave resonance. This is one of the 
reasons usually invoked to explain why these simplified models overestimate the scattering 
levels of Stimulated Brillouin Scattering (SBS). In this article we concentrate on SBS, which 
is the process by which the incident laser wave couples to an ion acoustic wave (lAW) to 
give rise to a scattered transverse wave. The generation of the harmonics due to the I AW 
fluid-type nonlinearity [1, 2, 3, 4, 5] is already known to be able to reduce signiflcantly 
the SBS reflectivity when compared with the results involving simply a linearized lAW. 
However, the previous fluid-type models for SBS in Refs. [1, 2, 3, 4], aimed at taking into 
account the JAW nonlinearity, were incomplete because they did not properly describe the 
flow modiflcation [6, 7] caused by the incident transverse wave momentum deposition. All 
the mentioned models [1, 2, 3, 4, 5] also ignored multi-dimensional effects. On the other 
hand, kinetic effects associated with particle trapping [8] give also rise to a nonlinear lAW 
frequency shift and therefore modify the SBS nonlinear behavior. 

In the present Letter, we reconsider the effect of the lAW nonlinearities on SBS by ac- 
counting properly for the flow modiflcation caused by SBS. We flrst derive approximate 
equations describing simultaneously the plasma hydrodynamics (i.e. the long wavelength 
density and flow proflles), SBS, and the harmonic generation of the excited lAW resulting 
from fluid-type nonlinearity. Our method consists in decomposing the fluid variables into 
long and short wavelength components, the latter corresponding to the SBS generated lAW 
and its harmonics[9]. Our new code, based on this harmonic decomposition method, makes 
it possible to describe plasmas of spatial sizes of the order of realistic laser produced plasmas 
(of mm-size, typically), because it does not resolve the lAW //m-scale. We then continued 
a step further by checking the capacity of our approach to account for kinetic effects ef- 
fects by implementing in the lAW propagator a nonlinear frequency shift modeling particle 
trapping. [8] 

The transverse electric field is described by £'(x, t) = e~*'^°* (^E^e^^°^ + i^.e"'*^"^^ + c.c. 
where £^+(x, t) and E_{-x,t) are the forward- and backward propagating hght field compo- 
nents, respectively, both enveloped in time and space with respect to the light frequency cuq 
and the wave number k^. This wave number is taken for a fixed reference plasma density 
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Neq which yields, using the critical electron density ric, — uJq{1 — Neq/n^/c^. For the 
plasma density n(x, t) and the velocity v(x, t) we use a decomposition separating the long- 
wavelength components A^o(x, t) and vo(x, t) and the short- wavelength components np(x, t) 
and Vp(x, t), with \p\ = 1, 2, 

n^No+ (nie''''''' + n2e'^'''' + + c.c.) , 
V = vo + (vie^'^^^ + Y2e'^'^ + .. + c.c.) , 

the first (p=0) representing the hydrodynamic evolution, and the terms with p >0, the 
fundamental ion acoustic wave, p =1, excited by SBS, and its harmonics, p >1. The 
reference wave number for the lAW is the wavenumber of backscattering, kg — 2ko, for 
which the ponderomotive force is proportional to oc E^E'^ exp{i2koz). 

We use the paraxial approximation to reduce the wave equation for the total electromag- 
netic field E to two "paraxial" equations for i5+(x, t) and £^_(x, t), 

CpUE+) = -iK/2nc) [niE_ +{No - N^)E+] , (1) 
CpUE_) = -i{uo/2n,) +(iVo - N,^)E_] , (2) 

with the paraxial operator £par(£^±) = [dt -\- c±dz + —i{c^ /2ujq)V'^^E±, where c+ and c_ 
stand for the group velocity of the forward/backward propagating light, respectively, with 
c+ = c^ko/cuo — — c_, and denotes the damping of the transversal waves. The right-hand- 
side (rhs) source terms in equations (1) and (2) account for (i) resonant 3-wave coupling due 
to SBS, with the fundamental ion sound wave, rii, and for (ii) refraction on long- wavelength 
density modifications, A^o ~ ^eq, causing e.g. self-focusing. In comparison with the full wave 
equation without decomposition into E±, this model allows a considerably coarser spatial 
resolution and thus much less numerical expense. 

For the long-wavelength hydrodynamic component we use the following set of equations, 
assuming isothermal conditions, and written in the conservative form on the left-hand side 
(Ihs): 

9t7Vo + V7VoVo = (9tnVw, (3) 
dt (iVoVo) + V (iVoVoVo) + c^VNo = (4) 

-NoclVUo + (atnv)iAw , 

where the rhs source terms, {dtn)iAw and {dtnv)iAw, describe the momentum transfer into 
the flow due to the lAW excitation by SBS, with {dtnv)^^^ = 2c5 {2usi — Vq • V) |nip/Ao, 
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and {dtn)ij^^ = — 2csV (|nip/A^o) ■ The ponderomotive force is given by VC/q = eoV(|-E+p + 
|£'_P)/ncTe. The equations describing the lAW driven by SBS, rii, and its harmonics, n/>i 
(using the convention n_; = n* for the complex conjugate) can be written, in the so-called 
weak coupling regime, as follows 



ni 



dt + J^si + m + (voz + Vgi)dz - i{cs/2lks)V 
-i{ksCs/2)No [8i,i (eoE+E*_/ncT) + 2Qi/N^] (5) 

with Qi = {l/2)Y.nhni-h for /i ^ and / ^ /i, where c, = [{ZTe + 3Ti)/Mi]i/2 is the lAW 
speed (with Z and Mi as the ion charge and mass), v^z the 2;-component of the flow vq; Vgi 
and uJi denote the group velocity and the "local" frequency of the /-th lAW harmonic, both 
accounting for the dispersion due to Debye shielding increasing with the harmonic order. 
They are given by Vgi — Cs(l + Z^A;^A|,)~^/^ and by 0Ji{z) — uJs{lksCs) + lksVoz{z) with the 
lAW frequency u;s{k) — kcs{l + /c^A|,)~-^/^. Equations (l)-(5) describe what we call the 
harmonic decomposition model. They form a closed system describing SBS in a temporally 
and spatially evolving plasma. They can be shown to conserve momentum [7] at the lowest 
order in l/{ksi\\) and in l/{ksi±y (with the inhomogeneity length £|| = \dzVo/vo\~^ and 
i± = |V±^;o/^;or')- 

In the following we emphasize the particular importance (i) of the SBS-induced flow 
modification, originating from the rhs term of Eq. (3) as well as of the term {dtnv)iAw 
on the rhs of Eq. (4), and (ii) of the lAW harmonic generation described by the coupling 
terms oc Y^nhTim-h in the rhs of Eq. (5). In order to stress the effect of each mechanism, 
we neglect for simplicity the lAW damping, (while being aware that the lAW damping 
coefficient is usually of the order of a few percent of the JAW frequency). Indeed, the 
SBS-induced flow modification due to momentum transfer, first pointed out by Rose in 
Ref.[6], cannot be ignored in the regime of absolute instabihty corresponding to weak lAW 
damping, because it is just in this regime that the stationary ID hmit of Eqs. (3) and (4) 
exhibits the most pronounced ffow modification. Namely, the generation of the backscattered 
light gives rise to a transfer of momentum to the bulk plasma in the spatial domain of 
SBS activity. This momentum transfer results in a decrease = vo,out — i'o,m < of 
the flow Vq in the direction of propagation of the laser, the net flow decrease being given 
by A^; ~ -2i?sBs(2eo|-B+| V-^eg7'e)(l - N^j2nc). Here, i?sBS denotes the SBS reflectivity 
corresponding to the considered SBS active region. 




FIG. 1: SBS reflectivity i?sBS versus time for the case of an undamped lAW with the parameters 
II = 2.5 ■ lO^^W/cm^ for Aq = 1.064/xm at = IkeV, iVo/^c =0.1 (taken at center), 2koXD =0.27, 
Lini ~ I6OA0. The soUd Hne is obtained from the decomposition code considering all terms, the 
dashed line from the decomposition code disregarding higher lAW harmonics, and the dash-dotted 
line from the "full" code. 



We have performed simulations on the basis of equations (l)-(5) and expanded the lAW 
up to its 3rd harmonic, resulting in a set of equations for rii, n2, and ris, with the rhs terms 
Qi — n2n\ + 713^2*, Q2 — n\-\- 2n^n\, and Qs = 3n2ni. We did not observe any significant 
changes when harmonics above the 3rd order were retained, while restricting to less than 3 
harmonics led to important differences. At this stage of our study we restricted ourselves to 
one-dimensional (ID) simulations in order to benchmark our harmonic decomposition code 
against a "complete" ID code which does not make the decomposition corresponding to 
Eqs. (l)-(3). This latter code solves Helmholtz's equation for the total electric field E{z,t) 
on the first hand, and the system of fluid equations for continuity and momentum, with the 
complete ponderomotive force, V\E{z,t)\^, as a source term, on the second hand. Here, in 
ID, the operator V reduces to the partial derivative e^dz. 

To ensure equivalent boundary and initial conditions we have considered a realistic case 
similar to an "exploding foil", where an initially heated plasma expands starting from an 
almost box-like density profile, with smooth shoulders, in the interval zi < z < Z2 along 
the laser axis. The plasma profile, with the initial plateau width Lini — I6OA0, successively 
undergoes rarefaction from each side, so that the velocity profile eventually tends to a 
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monotonous curve varying from negative to positive values with i>o = in the center. The 
simulation box is chosen in such a way that the rarefaction of the profile does not significantly 
change the boundary conditions for the light fields at the entrance {zent = < 2:1) and the 
rear side {z^ear > -22)- The total box size is Zj-^ar = 2000//i;o ^ 320Ao, where Aq = 271 /k^ 
denotes the laser wavelength. The boundary condition for the incident light at 2; = is a 
constant, £^+(0) = const, whereas the backscattered light is seeded with a noise source at 
the level = Zrear)!^) ~ 10~^|£'+(2; = 0)p and with a spectral bandwidth sufficiently 

larger than the lAW frequency, in order to cover all possible SBS resonances in the profile. 
In the density profile wings left and right of the central plateau (for times t < Lini/2cs), the 
plasma is strongly inhomogeneous in velocity and density so that SBS is inhibited by the 
strong flow gradient. 

We carried out our simulations in the absolute instability regime of SBS with undamped 
lAWs, both to examine the role of flow due to momentum transfer, and to benchmark the 
robustness of our decomposition code. Notice that in the case of completely undamped 
lAWs the SBS saturation level is, according to Refs. [10, 11], independent of the noise 
level. For the chosen electron density Ngq/ric = 0.1 and for the plasma length indicated 
above, the standard thrcc-wavc interaction model for undamped lAWs [10, 11] predicts a 
steep increase in the SBS reflectivities Rsbs as a function of the laser intensity, varying from 
RsBS « 1 for small laser intensities to i?sBS — 1 for normalized laser intensities above 
Cq = €o\Eo\'^/ncTe ~ 0.003, with £'0 denoting Eq — E+{z — 0). Our simulations comparing 
the decomposition code and the "complete" code show very good agreement, even for the 
extreme case shown in Fig. 1, corresponding to the plasma parameters mentioned above and 
to Og = 0.025 and 2koXD = 0.27 (corresponding to an electron temperature = IkeV and 
a laser intensity II = 2.5 ■ lO^^W/cm^ at Aq = 1.064/xm), for which the reflectivity would be 
99% in the absence of any JAW nonlinearity or flow modiflcations. For lower intensity values, 
the agreement is even more striking. This excellent agreement between the two codes gives 
us confldence in the robustness of the harmonic decomposition description. In the simulation 
presented in Figs. 1 and 2, the maximum amplitudes of the harmonics remained below the 
validity condition for harmonic expansion, namely \ni/Neg\ < Q^'^ilKXof for / = 1,2,3. 

It can be observed, in the spatial profiles shown in Fig. 2 for the backscattered inten- 
sity l-E-P, the fundamental lAW amplitude the flow Vq, and the plasma proflle A^o; 
that the lAW behavior and flow modiflcations are entirely connected with the existence of 
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FIG. 2: Spatial profiles, taken at 2koCst = 220, of the backscattered intensity (upper subplot), 
the fundamental lAW square amplitude |nip, the flow velocity vq, and the density Nq inside the 
shoulders of the exploding foil profile. 

"defects" in these spatial profiles and with a non-monotonous character in space (see also 
Ref.[4]). Namely, SBS develops in distinct spatial domains, interrupted by phase defects, 
which originate in the density profile shoulders corresponding to the low density plasma on 
the laser entrance side, and which then propagate into the profile plateau. Thus the SBS 
activity in each spatial domain appears to be uncorrelated, due to their different origin in 
the inhomogeneous velocity profile vo- This feature refiects in the structured nature of the 
backscattered light temporal spectrum, shown in Fig. 3 in which distinct peaks appear, and, 
consequently, in the temporally chaotic behavior of the reflectivity. 
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Our decomposition description makes it possible to discriminate the relative importance of 
the various effects contributing to the nonstationary behavior in the SBS reflectivity, as seen 
in Fig. 1. By suppressing parts of these effects in different runs, we have found that the most 
important effect is the excitation of the lAW harmonics: namely, retaining the harmonic 
excitation and neglecting the SBS-induced flow modiflcation lead to results that remain in 
reasonably good agreement with the exact model, whereas, retaining the flow modiflcation, 
but ignoring the harmonics leads to unphysically high levels of I AW amphtudes. It follows 
from these observations that a realistic modeling of SBS requires the proper description 
of the lAW harmonics. The ID simulations presented here would correspond to the SBS 
development in a long laser hot spot. We have recently carried out 2D simulations which 
conflrm the relevance in 2D of the scenario described above whenever the hot spot focus 
is not far (less than approximately one Rayleigh length) from the transition between the 
inhomogeneous and the homogeneous domain of the plasma density proflle (i. e. the shoulder 
of the expanding plasma in our case) . 

Increasing the laser intensity induces stronger lAW amplitudes at which ion and/or elec- 
tron kinetic effects take place. We have included phenomenologically weak ion kinetic 
effects in our decomposition model by adding a non-linear frequency shift of the form 
—ir}\ni/Neq\^/'^^ [8] in the propagator appearing in the Ihs of Eq. (5) describing the evo- 
lution of the lAW fundamental component and of its harmonics. Although this is subject of 
work in progress, let us mention that we have solved numerically the corresponding equation 
Eqs. (l)-(5), and we flnd that for a positive and sufficiently large 77 coefficient (~ 0.5 .. . 0.7), 
this shift can smooth out the effect induced by the harmonics and the flow, in a way such that 
(i) the "defects" are less pronounced and (ii) the SBS reflectivity diminishes, but without 
exhibiting a strong nonstationary behavior. In conclusion, we have shown that the SBS mod- 
eling presented here, based on a harmonic decomposition of the hydrodynamics variables, 
represent a promising way to describe laser plasma interaction in long scale-length plasmas. 
We have benchmarked our code based on the harmonic decomposition in the extreme limit 
of the absolute instability regime by neglecting the lAW damping. A novel picture of SBS 
arises in which an incoherent superposition of scattered light generated in distinct spatial 
domains in the velocity profile leads to a nonstationary character of the SBS reflectivity and 
to a significant reduction in the time averaged reflectivity. This harmonic decomposition 
description appears to be sufficiently robust and versatile to allow further sophistication by 
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FIG. 3: Spectrum of the backscattered light, corresponding to Fig. 1. The frequency is shifted 
with respect to the incident laser frequency loq. 



including additional mechanisms such as kinetic effects via an amplitude-dependent nonlin- 
ear frequency shift. We currently work on a generalization of the harmonic decomposition 
method in order to include the subharmonic lAW decay. 

The numerical simulation were carried out thanks to the access to the facilities of IDRIS 
at Orsay, France. The authors would like to acknowledge fruitful discussion with L. Divol, 
J. Myatt, C. Riconda, H. A. Rose, and W. Rozmus. 
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